Non-BCS behavior of optical properties across the cuprate phase diagram 
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The finite-frequency optical properties of tlie underdoped cuprates, in botli tlie normal and su- 
perconducting state, display features which go beyond a Fermi liquid and a BCS description. We 
provide an understanding of these properties within a simplified analytical model, which has been 
evolved out of the Hubbard model and ideas based on a resonating valence bond spin liquid. We 
find that: 1) in underdoped samples, the missing area integrals reveal a second energy scale due 
to the pseudogap, not present at optimum or overdoping; 2) the real part of the optical self-energy 
shows a large sharp peak, that emerges with the opening of the pseudogap which exists within the 
superconducting state and persists in the normal state; and 3) the amount of optical spectral weight 
which is transferred to the condensate is greatly reduced by the presence of the pseudogap as com- 
pared to the Fermi liquid case. These non-BCS features of the superconducting state are in good 
qualitative agreement with a body of experimental work on different cuprate systems and provide 
strong evidence from optical conductivity that they are all a manifestation of the pseudogap energy 
scale. 

PACS numbers: 74.72.-h,71.10.-w,78.20.Bh 



Understanding the approach to a Mott insulator from 
a Fermi Uquid state is an important question in cor- 
related electron physics. The cuprate superconductors 
provide a model system for investigating this issue. In- 
deed, they exhibit non-Fermi liquid and non-BCS behav- 
ior throughout an intermediate doping regime between 
the Mott insulator and the Fermi liquid. In the over- 
doped regime, the normal state appears to be rather con- 
ventional, though the superconducting ground state has 
d-wave symmetry. By contrast, the underdoped regime 
exhibits non-Fermi liquid features related to the pseu- 
dogap formation, which is thought to be responsible for 
many of its anomalous normal state properties. fl*! In ad- 
dition, in this regime, there is evidence that the pseudo- 
gap and superconducting correlations compete and both 
provide their characteristic imprints on superconducting 
properties, such as the two separate energy scales seen in 
Raman scattering. Q An important bulk probe of mate- 
rials is finite-frequency optical conductivity for which a 
considerable quantity of data has accumulated about the 
cuprates as a function of doping and in different systems. 
From this it is clear that the optical properties display 
anomalous behavior which is non-Fermi liquid and non- 
BCS-like. This behaviour has not been fully understood 
from a theoretical point of view nor has it been brought 
together to provide a conclusive statement about the ef- 
fect of a pseudogap energy scale in the superconducting 
state. In this commnication, we demonstrate that the 
non-BCS behavior can be understood as due to a second 
competing pseudogap energy scale within the supercon- 
ducting state. We provide three results in comparison 
with experiment which support this significant conclu- 



While several approaches to the modeling of the un- 
usual normal state of the cuprates have been proposed 
and examined, such as d-density wavesQ and phase in- 
coherent preformed pairs above TcQ, few have been able 
to address successfully the superconducting state or the 
evolution of properties with doping. Work on the Hub- 
bard (and related models), thought to be a good candi- 
date for describing the cuprate phase diagram, has been 
largely numerical but recently some progress has been 
made towards providing analytic approximations to these 
works which can facilitate ease of calculation and provide 
a deeper understanding of the underlying physics. In or- 
der to examine the issue of the optical conductivity and 
its non-BCS behavior, we have chosen to adopt an analyt- 
ical approach proposed by Yang, Rice and ZhangQ which 
is evolved out of work that has its basis in the Hubbard 
modelQ and the ideas of a resonant valence bond (RVB) 
spin liquid, proposed by Anderson^, 0]. The coherent 
part of the electronic Green's function in this theory en- 
codes the antiferromagnetic Brillouin zone with pseudo- 
gap formation, as the Mott transition is approached. We 
use this formalism to calculate the in-plane optical con- 
ductivity a{ijj) in order to provide an understanding of 
the anomalous superconducting behavior seen in quanti- 
ties such as the partial optical sum and the optical self- 
energy. 

In the analytical model used here, the coherent part of 
the electronic Green's function takes the form[5|] 
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FIG. 1: (color online) (a) Maximum value of the pseudo- 
gap and superconducting gap in units of to versus doping (x) 
dependence:^] A°g{x)/2 = 0.3(1 - x/0.2) and A°^{x)/2 = 
0.07[1 — 82. 6(1 — 0.2)'^]. (b) Luttinger area in the pseudogap 
state for four dopings indicated with arrows in (a), kx and ky 
are in units oi n/a, with a the lattice spacing. G(k, 0) > for 
the shaded region, the area of which satisfies the Luttinger 
sum rule. 



with self-energy 
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The pseudogap self-energy is given by I]pg(k, w) = 
Apgp/[a; -f- eo(k)]. In these formulas gt is the relative 
weight of the coherent part of the Green's function, which 
we set equal to one. There is also an incoherent piece not 
included here. The superconducting gap Asc and pseu- 
dogap Apg are both assumed to have d-wave symmetry, 
i.e. A(k) — A°(a;)[cos(fc2:a)—cos(fcj^a)]/2, within the two- 
dimensional square Cu02 Brillouin zone. The doping (x) 
dependence of their amplitude A°(x) is given in Ref. [5| 
for a particular case and is reproduced here in Fig. [T^. 
Finally, e(k) is the electronic dispersion and eo(k) the 
first- nearest- neighbor-only version, with eo(k) = defin- 
ing the antiferromagnetic Brillouin zone. In Fig. [T]3, we 
show the Luttinger surfaces of the theory for the param- 
eters of Ref. The shaded blue areas correspond to 
occupied states containing an area proportional to 1 — a: 
and the white pockets within the first antiferromagnetic 
Brillouin zone enclose an area proportional to a;. To com- 
pute the in-plane optical conductivity one needs to know, 
in addition to the Green's function, the Gorkov anoma- 
lous amplitude denoted by F'^ (k, lo) and given by 
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From these, the conductivity (t{T,lu) follows from the 
Kubo formula 

^m-) = ^E^k r f^[/H-/(c.+.)] 

X [A{k, oj)A{k, uj + i^) + B{k, uj)B{k, oj + ly)^) 

where f{uj) is the Fermi-Dirac distribution function at 
temperature T, Wk is the electron velocity, and the spec- 
tral functions are given as A(k, ui) = — 2ImG(k, ui + iO~^) 




FIG. 2: (color online) (a) and (b): Optical conductivity curves 
with (solid red curve) and without (dashed black curve) a 
pseudogap. (c) and (d): Plots of A(k,tj = 0) to show the 
Fermi arcs, with kx and ky labels normalized as indicated in 

Fig.m 



and i3(k, Lj) = — 2Imi^^^(k, + zO+). We present results 
at T = for the generic band structure and parame- 
ters suggested for cuprate systems @. In Fig. [21 we show 
ai{uj)/A as a function of photon energy uj/to, where to 
sets the energy scale of the first nearest neighbor hopping 
parameter. The normalization A has been defined as the 
area under the conductivity curve in the normal state 
for X — 0.2 (optimally doped). Fig. ^jp. is for x — 0.05 
(highly underdoped) and Fig. [2b is for a: = 0.1, just in- 
side the inner edge of the superconducting dome as shown 
in the phase diagram of Fig. [1^. In Fig. [2fc and d, we 
plot the corresponding color maps of ^(k, = 0) which 
provide information about the Fermi arcs as a function 
of doping. In Fig. ^ and b, the dashed black curves 
contain no pseudogap and are similar to a Drude form, 
while the solid red curves contain a finite Apg. In all our 
calculations we have replaced the Dirac delta functions 
in the spectral densities A and B with Lorentzians of 
width O.Olio. Twice this value gives the corresponding 
transport width 1/t. The solid red curves, in contrast to 
the dashed curves, deviate strongly from Drude behavior. 
The optical spectral weight in the Drude-like distribution 
centered at = is greatly reduced, and some, but not 
all is transferred to the interband transition region. The 
onset of the region to which this weight is transferred, is 
close to for x = 0.1, but falls considerably short of 
Apg for X = 0.05. The reason for this is that the gap value 
on the Luttinger pockets is much smaller than A^^ (see 
Fig-[2fc)- For both dopings, the rapid drop denoting the 
end of the interband region is located at approximately 
2Apg. These onset energies can be traced to peaks in the 
electronic density of states. Not all the optical spectral 
weight lost in the Drude region is recovered in the in- 
terband region at higher energies. We are dealing here 
only with the coherent part of the Green's function. In a 
more complete theory some of the missing weight could 
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FIG. 3: (color online) (a) and (b): Optical conductivity curves 
with (solid red curve) and without (dashed black curve) a 
superconducting gap. The pseudogap is also present in (a) 
as can be seen in the phase diagram in Fig. [1] (c) and (d): 
A(k,tij = 0) showing the Fermi arcs (labels same as Fig. 



FIG. 4: (color online) In (a) and (b) the shaded region rep- 
resents the missing area, that is the spectral weight lost to 
the condensate for x — 0.14. The normal state in panel (a) 
does not include a pseudogap, whereas in (b) it does. Panel 
(c) shows the integral of the missing area S(w) normalized 
to ^(co). The inset in (c) shows the percentage of spectral 
weight lost to the condensate as a function of doping for a 
normal state without (blue circles) and with (red crosses) a 
pseudogap. 



be transferred from coherent to incoherent part. Similar 
to our result, a two-component optical conductivity has 
also been found in numerical work[9] on the t — J model. 

In Fig. [3] we compare the normal state (dashed black 
curve) with the corresponding superconducting state 
(solid red curve) for (a) x = 0.14 and (b) x = 0.2. In 
this latter case, Apg = so that the corresponding Lut- 
tinger surface reduces to the usual large Fermi surface 
of Fermi Hquid theory (Fig. [3Ji), and the real part of 
the conductivity, cri(w), in the superconducting state be- 
haves in a conventional way. For an s-wave supercon- 
ductor there would be a gap of 2A^^ in ai{ui) with all 
optical spectral weight below this energy going into the 
condensate and at higher energies, the normal state con- 
ductivity is recovered. However, here the gap has d-wave 
symmetry, and some weight remains in the Drude re- 
gion below 2AO^.[i3] Fig.lia is more unusual because the 
pseudogap, which exists in the normal state, is responsi- 
ble for depleting the Drude peak and transferring weight 
to the interband region. Adding superconductivity fur- 
ther depletes the Drude region below 2A^^, although it 
does not eliminate it completely. In addition, and this 
goes beyond the BCS model, it shifts the onset of the 
interband transitions to higher energies. The differences 
between the cases with (frame (a)) and without (frame 
(b)) pseudogap are further emphasized in Fig. [4l where 
we show the real part of the conductivity for x = 0.14. 
The residual scattering rate is 2r/ = 0.02to in both cases. 
The shaded regions show the missing area under ai{uj) 
due to the superconducting transition. Frame (a) is con- 
ventional d-wave BCS behavior but frame (b) is anoma- 
lous. Comparing between (a) and (b), below w/to ~ 0.2, 
erf (oj) is not affected much by pseudogap formation be- 
cause the remaining absorption (Drude) only involves the 
nodal direction. However, in this energy region (uj) is 
much larger in (a) than it is in (b) because in this lat- 



ter case the antinodal direction is gapped out with its 
spectral weight transferred to the interband transitions. 
Note that this interband piece contains about 20% of the 
condensate. This is further elaborated upon in (c) where 
we show results for the missing area partial optical sum 
S{uj) = jQ+[cr^{j^) — ai{v)]dh'. The dashed blue curve is 
without, while the solid red curve is with a pseudogap. 
Pseudogap and superconducting energy gap scales are in- 
dicated by arrows. The dashed blue curve, normalized to 
S'(oo), rises rapidly to its saturation value on the scale 
of the superconducting gap for the particular residual 
scattering rate used (a quantity also involved in deter- 
mining the rate of initial rise). For the solid red curve, 
this first rise is also clearly visible but in addition there 
is a second region of slower rise on an energy scale set by 
the pseudogap which involves about 20% of the conden- 
sate. This behavior has been qualitatively verified in the 
experimental work of Homes et al. 11] for YBa2Cu306.95 
(no pseudogap}and YBa2Cu3O6.60 (with pseudogap) [see 
Fig. 7 of Ref. ll|]. A quantative comparison is not possi- 
ble because in the analysis of the experimental data, the 
normal state at low temperature is not available and its 
value just above Tc is used instead. Still the pseudogap 
scale is seen only in the underdoped case. 

In the inset of Fig.jH;, we show the fraction of the total 
available optical spectral weight that goes into the con- 
densate as a function of doping. For an ordinary s-wave 
BCS superconductor this amount would be 100% in the 
clean limit, l/r ^ A°^, as it would also be for a d-wave 
gap. The entire optical spectral weight condenses and 
the zero temperature penetration depth depends only on 
normal state parameters and not on the value of the gap. 
For finite 1 /r this no longer holds and the amount in the 
condensate is reduced with the reduction depending on 
the relative size of l/r to A°^. However, when a pseu- 
dogap is introduced, the condensate fraction is further 
reduced. This can be traced to the transfer of spectral 
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FIG. 5: (color online) (a) Results for —2Yf-^ {to) vs lo from 
the theory (only x = 0.2 has no pseudogap). (b) The data 
of Hwang et al.^ for the A{-2E°''(a;)} vs lo in cm"! for 
underdoped Bi2212 with Tc — 67K (an incoherent background 
has been subtracted out). The black dashed-dotted curve is 
the theory for to ~ 450 meV. Both (a) and (b) show the 
superconducting state. 



weight from the Drude to the region above A^g where it is 
not as susceptible to condensation. Note in particular the 
points a,t X — 0.14. The open circle without pseudogap 
is close to 70% while the opening of a pseudogap reduces 
this value to roughly 40%, less than half the electrons 
participate in the condensation. This is an important re- 
sult associated with the non-BCS behavior of this theory 
and for which there is experimental evidence. ^12] 

Finally in Fig. [5l we present our results for the real 
part of the optical self-energy T,°p{uj). It has become 
common in correlated systems to describe the conductiv- 
ity a{uj) in terms of a generahzed Drude form[l^ with 
cr(r,w) = {iujl/4:Tr)/[uj - 2I]°P(r,w)], where ujp is the 
plasma frequency. For the identification of the pseudo- 
gap in in-plane optics it was recently found that S^''(T, lu) 
is a more useful quantity Q than is the conductivity it- 
self, and it shows a prominent "hat" -like structure seen 
above a large inelastic background. Prominent peaks in 
Y^'i'{T,ll}) can arise from either pseudogap or supercon- 
ducting gap formation 15, [l^ when the inelastic scat- 
tering, described by the incoherent part of the Green's 
function, involves sharp boson structure on the same en- 
ergy scale. In the present work such structures arise even 
though inelastic scattering is not accounted for. While 
it is not clear exactly how to separate coherent and in- 
coherent contributions, Hwang et al.[l6j have removed 
a background from their data presented in their Fig. 19 
for the BiaSrzCaCuzOs+s (Bi2212) series. The results 
of our calculations, which are shown in Fig. [5^ for var- 
ious values of doping x, are best compared with these 
data for one of the most underdoped samples studied 
with Tc = 67K, reproduced in Fig. [5)3 (sohd blue curve). 
Pseudogap effects are large in this case and dominate 
the peak formation, which is also hardly affected by su- 
perconductivity in both experiment and theory. Tak- 
ing to = 450 meV and x — 0.15, the black dot-dashed 
curve is reproduced in Fig. [SJs. The agreement with data 
is reasonable particularly as we have not used param- 



eters specific to Bi2212 but instead stayed with those 
for Ca2-a;Naa;Cu02Cl2 given in Ref. [Hj. The subtrac- 
tion of the incoherent contribution to S]°^(r, lu) may also 
be responsible for some of the discrepancies seen beyond 
~ 900 cm~^. As stated the peaks of Fig. [5^ are greatly 
reduced with increasing doping, a trend confirmed in the 
normal state experimental data. 

In summary, we have been able to explain three ma- 
jor anomalous experimental observations from finite- 
frequency optical conductivity measurements on several 
different cuprate systems where non-BCS behavior is ex- 
hibited as a function of doping. Such non-BCS behavior 
is understood in our work to result from a pseudogap- 
type energy scale present in the normal state which also 
manifests itself within the superconducting state. Com- 
parison of our theoretical results with optical data in the 
Bi2212 series reveals good qualitative agreement for the 
evolution with doping of both position in energy and size 
of the normal state "hat" -like structures seen in the real 
part of the optical self-energy. In addition, the partial 
optical sum reveals that a second energy scale related to 
the pseudogap exists in underdoped samples and is not 
present in optimally or overdoped samples, as seen in 
the data of Homes et al[llj| for YBa2Cu307_5. Another 
significant result of our analysis is that the opening of 
a pseudogap greatly reduces the percentage of carriers 
that condense in the superconducting state as compared 
to BCS. We stress that optical conductivity is a bulk 
probe applicable to many cuprate systems, as opposed 
to surface probes such as scanning tunneling and angle- 
resolved photoemission spectroscopies which are mainly 
limited to the Bi2212 system. Consequently our analysis 
provides strong support that the observed non-BCS be- 
haviors seen in several cuprate systems and their doping 
dependence can be described via a theory, such as we 
have used here, based in the Hubbard model and ideas 
of a RVB spin liquid. 
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